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> ■ 1 INTRODUCTION 



ABSTRACT 

We present the results of a comparative period search on different timescales and 
modeling of the ZZ Ceti (DAV) star GD 154. We determined six frequencies as normal 
modes and four rotational doublets around the ones having the largest amplitude. Two 
normal modes at 807.62 and 861.56 /iHz have never been reported before. A rigorous 
test revealed remarkable intrinsic amplitude variability of frequencies at 839.14 and 
861.56 /iHz over a 50-day timescale. In addition, the multi-mode pulsation changed 
to monoperiodic pulsation with an 843.15 /xHz dominant frequency at the end of the 
observing run. The 2.76 /^Hz average rotational split detected led to a determination 
of a 2.1-day rotational period for GD 154. We searched for model solutions with 
effective temperatures and \ogg close to the spcctroscopically determined ones. The 
best-fitting models resulting from the grid search have Mh between 6.3 * 10^^ and 
6.3* 10~^M*, which means thicker hydrogen layer than the previous studies suggested. 
Our investigations show that mode trapping does not necessarily operate in all of the 
observed modes and the best candidate for a trapped mode is at 2484 /^Hz. 
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It has been well established that g-mode white dwarf pul- 
sators show diversity in light variation from the simple si- 
nusoidal to the mightily complicated ones. The latter vari- 
ation accompanied by nonlinear features due to harmon- 
ics and combination freq uencies (see e.g. the review of 
iFontaine fc Brassardl bOOSTl. 

It has also been pointed out that pulsational amplitudes 
and/or the frequency content of light curves might vary 
on short timescales. These phenomena are quite common 
amongst the cooler DAVs and DBVs, and also the PNNV 



stars (|Handler..2003. 1 . The amplitude variability can be help- 
ful to asteroseismology as different normal modes become 
detectable over time and therefore the n umber of known 
modes increase (see e.g. the case of G29-38. lKleinman et al.l 
119981 ). However, amplitude variations on timescales of weeks 
and months can make the accurate determination of indi- 
vidual modes difficult or even impossible using the stan- 
dard Fourier deconvolution technique. These variations oc- 
cur sometimes as sudden effects as in the case of the 'ev- 
erchanging' GD 358, the most spectacular representative of 
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the phenomenon (JKepler et al.ll2003l : IProvencal et al.ll2009l ) . 
In spite of intense observational efforts, the origin of these 
phenomena still remains unknown. A possible explanation, 
which can be tested relatively easily, is mode beating due to 
different unresolved pulsational modes (such as rotational 
splitting) . 



The DAVs lying near the red edge of the ZZ Ceti in- 
stability strip (GD 154 is an example) are characterized by 
long periods and complex pulsational behaviour. Changes 
in pulsation from a single mode status (connected by har- 
monics and not trivially by subharmonics) to a multi-mode 
status have been seen for GD 154 during the previous obser- 
vations. Since the seemingly new appearance of a different 
normal mode can be caused only by the amplitude varia- 
tion of the ever excited normal modes, we decided to follow 
the pulsational behaviour for an observational season. We 
present our comparative analyses on different timescales. It 
was the increased number of pulsational periods determined 
by our dataset that allowed us to invstigate the star from an 
asteroseismological point of view. We also discuss the related 
results in the paper. 
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Table 1. Log of observations of GD 154. Five subsets were created 
based on the closely spaced nights. 



Run 


Subset 


UT date 


Start time 


Points 


Length 


No. 


No. 


(2006) 


(BJD-2 450 000) 




(h) 


01 


1 


Feb 03 


3769.520 


372 


4.04 


02 


1 


Feb 05 


3771.575 


335 


3.28 


03 


1 


Feb 07 


3773.519 


340 


4.06 


04 


2 


Mar 02 


3797.430 


585 


5.66 


05 


2 


Mar 06 


3801.481 


484 


4.72 


06 


2 


Mar 07 


3802.285 


787 


9.11 


07 


2 


Mar 08 


3803.305 


684 


6.60 


08 


3 


Mar 31 


3826.364 


659 


6.60 


09 


3 


Apr 02 


3828.288 


291 


5.03 


10 


3 


Apr 04 


3830.357 


451 


4.42 


11 


4 


Apr 21 


3847.385 


493 


5.45 


12 


4 


Apr 22 


3848.321 


709 


6.88 


13 


4 


Apr 23 


3849.316 


549 


6.85 


14 


4 


Apr 24 


3850.380 


468 


5.20 


15 


4 


Apr 25 


3851.332 


541 


6.44 


16 


5 


Jul 13 


3930.340 


65 


0.71 


17 


5 


Jul 17 


3934.335 


92 


1.06 


18 


5 


Jul 18 


3935.336 


161 


1.83 


19 


5 


Jul 19 


3936.342 


150 


1.65 


Total 








8216 


90.19 





Figure 1. The CCD field with the variable and the comparison 
stars. The marked stars CI, C2, C3 were used to construct a 
reference light curve. 



Figure 2. Differential light curve of CI to the average of the three 
reference stars is shown (CI, C2, C3), obtained on JD 2 453 803. 
The standard deviation is 0.003 mag. 

The original images were reduced using standard lRAl{j 
routines. Bias, dark and flat corrections were applied on the 
frames. We performed aperture photometry using the iraf 
DAOPHOT package, setting the aperture size to two times the 
average FWHM on the given night, and applied differential 
photometry on the stars. All the time data were calculated 
to Baricentric Julian Date (BJD). 

We obtained BVRI photometry and checked the B — V 
indices of the potential comparison stars on the frames, but 
we found that they were rather different than the value of 
GD 154 {B — V = 0.18). As there was no candidate for a 
comparison star with similar colour to our target, we tested 
all the stars on our CCD field by computing the average 
magnitudes to construct a reference hght curve. We found 
that the most constant signal with the lowest standard de- 
viation value could be constructed by averaging the data of 
the three brightest stars around the variable (see Fig. [T|. 
A part of a differential light curve of the reference system 
and a check star's (CI) averaged data is presented in Fig. (2] 
We applied the widely used polynomial fitting method to 
remove the effect of the atmospheric extinction. The low- 
frequency part of the Fourier Transform (FT) is biased by 
this filtering (up to 20-30 cd~^), for which reason we did not 
focus on low frequencies in our Fourier analysis. 



3 FREQUENCIES IN GD 154 

The hght curves of GD 154 during 2006 can be seen in Fig.O 
The different minimum to maximum amplitude of the cycles 
suggest an overall multi-mode behaviour, except for the reg- 
ular cycles of the last nights. To get an overall view of the 
pulsation modes of GD 154, we present here the modes de- 
termined in different epochs. 



2 OBSERVATIONS AND DATA REDUCTION 

We used the 1-m Ritchey-Chretien-Coude telescope at 
Piszkesteto mountain station of Konkoly Observatory for 
data collection on GD 154. The observations were made with 
a Princeton Instruments Vers Array :1300B back- illuminated 
CCD camera without any filter. Our observing season 
spanned six months in 2006. Altogether, 90 hours of pho- 
tometric data were collected on 19 nights. 30 s integration 
times were used. The longest contiimous single-inght hght 
curve was more than 9-hour long. The shortest one that we 
found usefijl to determine the frequency content covers two 
full cycles. Details of the observations are provided in Ta- 
ble [1] 



3.1 Frequencies in the previous data 

The hght varia t ion o f GD 154 was discovered in 1977 by 
I Robinson et al.l ([1973). The star showed nonlinear monope- 
riodic pulsation with F = 843 piHz main frequency, harmon- 
ics (2F, 3F, 4F, 5F) and 'intermediate frequencies' (1.52F, 
2.53F, 3.53F). The pulsation was regular except on the last 
night of the observations, when the hght curve changed dra- 
matically as the 1.52F mode became dominant. The authors 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 4. Fourier spectrum of the wliole dataset. Tiie window 
function is given in the insert. 



pulsation. Three independent modes (/i = 842.8 pHz, /2 = 
918.6 /iHz, /a = 2484.1 /xHz) and their triplet components 
were found in the dataset. All the other peaks in the power 
spectra were explained as linear combinations and harmon- 
ics. Non-linear behaviour was clearly confirmed but no dras- 
tic period change was observed. Surprisingly, they did not 
find peaks around half-integer frequencies ijPfeiffer et al.l 
Il996l ). 

In 2004, an analysis of a two-site observational cam- 
paign (JHiirkal et al.ll2005l ) reported changes in the frequency 
and amplitude content. Some additional frequencies ap- 
peared, however none was near the subharmonic. They were 
interpreted as new excited modes at 786.5 yRz, 885.4 yRz 
and 1677.7 /xHz. However, the authors did not find the fre- 
quencies at 842.8 /iHz and 918.6 /iHz that were observed be- 
fore. 

We can summarize that although four modes (786.5, 
842.8, 885.4 and 918.6 /iHz) were reported in a 132 ^Hz 
range, not more than two were present during a given observ- 
ing run. None of them was present in every d ataset, although 
one o f them (Ri843 /iHz) was presen t in 1977 (|Robmson et al.l 



119781 ) and in 1991 (WET campaign. [Pfeiffer et al.ll 19961 ) with 
nearly the same value. The different solutions raise the pos- 
sibility that (1) different modes are excited from time to 
time or (2) their amplitude is changing to below or above 
the detection limit, or (3) the actual frequency content in- 
terfering with the alias patterns of a given data distribution 
results in a seemingly changeable frequency content. 



Figure 3. Normalized differential light curves of the 19-night run 
on GD 154. 



suggested that the energy exchanged between the two modes 
due to weak non-linearity. 

The existence of the near half-integer frequencies 
in the power spectra is not unique. Similar pulsa- 
tional behaviour was observed in PG 1351-1-489 (DBV, 
Goupil. Auvergne fc Bagiiiil Il988l ) and in G191-16 (DAV, 
Vauc lair ct al. 1989|). According to a possible explanation 
an independent g-mode can be excited near 1.5F, as we see 
in the interpretation of the data obtai ned on BPM 31594 
l|0'Donoghue. Warner fc Croppeij[l993 ). 

The Fourier spectral features and the sudden change 
in the dominant period made GD 154 interesting enough 
to be the target of the Whole Earth Telescope (WET, 
iNather et al.l Il990l ) campaign in 1991. The 12-day quasi- 
continuous observations and the follow-up campaign orga- 
nized a month later resulted in another interpretation of the 



3.2 Analyses of the new data 

A standard frequency analysis was performed with 
the Multi-Frequency Analyzer (MuFr An) time-string tool 
l|Kollathlll99"ol : lcsubrv fc Kollathll2004l ). The software pack- 
age is a collection of methods for period determination, sine- 
wave fitting for observational data and graphic routines for 
visualization of the results. In each step, a fine tuning of the 
actual multi-frequency solution was carried out to avoid the 
disadvantage of prewhitening. At the same time, both the 
amplitudes and the fit were checked to avoid any misidentifi- 
cation of a peak as a real frequency. Investigations for signif- 
i cance and errors wer e obtained with the Period04 package 
dLenz fc Bregerll2005l ). 

The Fourier spectrum of our whole dataset (over 167 
days) revealed the complexity of the resolution based on a 
set of single site observations (Fig. |4|. 

The overall appearance, a pyramid-like feature, suggests 
that more frequencies are present in the 760 - 960 /xHz fre- 
quency range. The fine structure is more complicated than 
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Table 2 . The frequency content of the whole data and two subsets 



Mode 


Whole data 


2nd subset 


4th subset 




Freq. Ampl. 


Freq. Ampl. 


Freq. Ampl. 


No. 


/iHz mmag 


/jHz mmag 


fiRz mmag 



1 


839.14 


9.16 


838.92 


16.19 


838.53 


6.00 


2 


843.15 


9.46 


- 


- 


843.50 


4.07 


3 


844.65 


7.22 


845.17 


3.82 


- 


- 


4 


861.56 


7.31 


- 


- 


861.52 


13.04 


5 


864.27 


4.23 


- 


- 


863.98 


4.39 


6 


857.84 


4.37 


- 


- 


- 


- 


7 


807.62 


3.99 


- 


- 


- 


- 


8 


802.94 


3.97 


- 


- 


802.93 


5.81 


9 


809.20 


3.76 


809.77 


8.33 


810.17 


3.18 


10 


918.72 


4.78 


918.95 


5.93 


917.80 


5.26 


11 


920.41 


4.11 


- 


- 


920.24 


7.10 


12 


883.56 


4.32 


888.06 


2.80 


883.48 


4.44 



0.015 
0.01 


1 1 1 2nd subset " 


■3 0.005 

CO 

M. 


i^,,,,,v;.ytyri^ 1 '" "I' 1 mkA^nrnhfyh^^^ 
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Figure 5. Comparison of the frequency content of the 2nd and 
4th subsets. The frequency content of GD 154 and especially the 
amplitudes changed from 3803 to 3847 BJD. 



the spectral window (insert in Fig. |4| and suggests that the 
amplitude of the modes may have changed during the whole 
timebase. 

We carried out comparative analyses on different 
timescales (weeks and a whole observational season) that re- 
sulted in different spectral window patterns. On a medium- 
size timescale (regularity of telescope allocation) the fre- 
quency resolution is acceptable, and the amplitude variation 
is still not masked. The baseline of the whole dataset gives 
a better resolution, but the amplitudes are not correct when 
amplitude variation is going on in a shorter timescale. 

Concerning the 760 - 960 /iHz frequency region, twelve 
frequencies were identified in the whole dataset, and are 
given in Table [2] (col. 2 and 3). The 0.035 /.iHz Rayleigh fre- 
quency of the whole dataset guaranteed that the frequencies 
are properly resolved. According to the frequency differences 
three triplets, a doublet and a single frequency were found. 
We confirmed the frequenci es found prev i ously, except the 
one at 786.45 ^Hz given bv iHiirkal et al.l l|2005l ). Our new 
discoveries are two modes at 807.62 and 861.56 /xHz values. 
The increased number of normal modes puts more constrains 
for modelling. 

Checking the time dependence of the frequency content, 
five subsets were created, but only the longest and the ob- 
servationally most populated subsets' Fourier spectra (2nd 
and 4th) are presented here (Fig. [S]). In the 2nd subset (four 
nights between 3797 and 3803 BJD) 26.1 hours of observa- 
tions were collected over 5.9 days. In the 4th subset (five 
nights between 3847 and 3851 BJD) more measurements 
were obtained (30.8 hours) but over a shorter timebase (3.9 
days). The corresponding Rayleigh frequencies on the time- 
base of the subsets are IpHz and 1.57 /iHz. The alias pat- 
tern of the 4th subset is much cleaner due to the consecutive 
nights. 

At first sight the Fourier spectrum of the 2nd and 4th 
subsets (upper and lower panels of Fig. [5]) suggests that the 
frequency content of GD 154 and especially the amplitudes 
changed from 3803 to 3847 BJD. 

We clearly recognize distinct peaks at 839 /iHz and 
918 ^Hz in the 2nd subset and at 861 ^Hz and 802 fiHz in the 
4th subset. The accepted frequency content of the two sub- 
sets as a result of comparative analyses are given in Table [2] 



We found a doublet around the dominant mode in the 2nd 
subset and four doublets in the 4th subset, for a total of 5 
and 9 frequencies, respectively. The triplets are not resolved 
on the subsets. Although the frequency values obtained on 
the different datasets agreed within 1 /iHz in most cases, the 
amplitudes are remarkably different at different epochs. 



3.2.1 Results on data simulation 

Two points were checked in synthetic data. The fre- 
quency spacings of 839.14 /xHz and 883.56 /tHz frequencies 
to the newly found mode at 861.56 /iHz are 22.42 /iHz and 
22.0 /iHz, about twice the daily alias pattern (11.57 /tHz), 
that could be confusing in the case of a single site observa- 
tion. The reason of the apparently variable amplitudes of the 
839.14 /iHz and 861.56 /tHz frequencies in the two subsets 
were also checked, whether they are interaction of the unre- 
solved triplets or they represent intrinsic amplitude variation 
of normal modes. 

Synthetic data were generated with the different combi- 
nation of the possible normal modes, doublets and triplets, 
with constant amplitude for the time series of the whole sea- 
son and the two subsets. The analyses of the synthetic data 
were conclusive. 

The five possible normal modes can be obtained with 
high precision not only for the whole dataset but for the 
subsets, too. We confirmed that the new frequency at 
861.56 /iHz is a real normal mode, it cannot appear by the 
daily alias interaction. The interaction of doublets does not 
cause large amplitude changes in the subsets. The interac- 
tion of the three triplets on shorter timebase can cause as 
large amplitude increase as the one we found in the analy- 
ses of our real time series. However, the amplitude ratio of 
the 838.92 /iHz and 861.52 /tHz frequencies did not change in 
the synthetic data. The amplitude of 838.92 /iHz frequency 
was always higher in both subsets than the amplitude of 
the 861.56 /tHz frequency. We confirmed that both modes 
reveal real amplitude variations over a 50-day timescale. As 
we require a stable frequency solution for our single site 
measurements, we present a solution that includes only the 
doublets. 
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folded with 842.93 mikroHz 



Table 3. The frequency content of GD 154 at 2006. Four doublets 
and two independent modes were found. The modes at 861.56 and 
807.62 ^Hz have never been reported before. 
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Figure 6. The light curves of the 5th subset (BJD 3930 - 3936) 
folded by the dominant mode of the subset (upper panel) and 
by the dominant mode of the whole dataset (lower panel). The 
regular arrangement of the upper panel suggests a monoperiodic 
pulsational phase of GD 154 during the given interval. 



3.2.2 Harmonics and subharmonics 

We performed frequency analyses in the region of harmon- 
ics and subharmonics, too. An independent frequency at 
2484.14 ^Hz was definitely found on those nights when the 
star showed multimode pulsation. 

We found groups of peaks around 1652.78, 1711.81, 
1752.31 and 1836.81 /xHz values but they do not correspond 
exactly to twice the value of the mother frequencies and it is 
hard to explain them as linear combinations. However, the 
second harmonic of the dominant mode (2517.48 pHz, the 
mother frequency is 839.14 /iHz), near to the high frequency 
independent mode (2484.14 //Hz), was clearly recognized. 

Linear combinations of the dominant mode and its dou- 
blet with the high frequency independent mode were de- 
tected at 3253.47 and 4155.09 /iHz. Although these values 
are near the 4F and 5F values, they agree with the linear 
combination explanation much better. 

On the last, but unfortunately the shortest nights only 
a single mode at 842.93 //Hz seems to be excited, instead 
of a multi-mode behaviour. It corresponds to the rotational 
component (at 843.15 /iHz) that was dominant in the dis- 
covery and the WET runs. In Fig. |6] we present the light 
curves of the 5th subset (BJD 3930 - 3936) folded by the 
dominant mode of the subset, 842.93 fiHz (upper panel) and 
by the actual value of the dominant mode of the whole 
dataset, 838.94 /iHz (lower panel). The narrow spread of 
the measurements in the upper panel proves that the star 
has a single mode dominant and this dominant mode is 
842.93 AtHz in the BJD 3930 - 3936 interval. Unfortunately, 
we could not follow in details this status of the pulsation 
that seems to appear from time to time between the multi- 
mode states. On 20 and 21 April, 2007 GD 154 showed a pure 



Freq. 


Period 


Ampl. 


Phase 


(/.Hz) 


(s) 


(mmag) 


(degree) 


839.14 


1191.7 


8.38 


314.4 


843.15 


1186.0 


6.90 


1.9 


861.56 


1160.7 


5.76 


276.5 


864.55 


1156.7 


4.91 


152.8 


918.70 


1088.5 


4.44 


215.6 


921.61 


1085.1 


3.88 


94.6 


807.62 


1238.2 


4.64 


218.9 


803.74 


1244.2 


4.38 


144.0 


883.56 


1131.8 


4.12 


337.4 


2484.14 


402.6 


3.50 


82.6 



monoperiodic pulsation state (Montgomery, private commu- 
nication) with harmonics and subharmonics. The convec- 
tive response timescale of the DAV star EC14012-1446 was 
compared partly to that val ue of GD 154 ob tained from 
the monoperiodic light curves (jProvencal et al.ll20 12l. Weak 
signs of the subharmonics (at 1195.60 and 1291.67 //Hz) ap- 
peared mostly on our last short nights connected to the sin- 
gle mode pulsation. There were nights when the sign of the 
2.5 times (2104.17 /iHz) or 3.5 times (3003.47 /iHz) values of 
the mother frequency could also be found. 

The Fourier parameters of the final accepted frequency 
content of GD 154 are given in Table |3] 



3.2.3 Characteristic features m the light curves 

A fit, generated by the Fourier parameters given in Table |31 
describes the general features of the light curve of GD 154. 
However, there are intervals containing some (4-6) cycles 
where the measurements represent much higher amplitudes 
for the cycles than the fit. Two special features were iso- 
lated. In panel a of Fig. [T] the envelope of the measure- 
ments resembles a Gaussian profile with steep increase be- 
fore the maximum and a steep decrease after it. The fit by 
our mode decomposition gives a much lower regular ampli- 
tude of the cycles. Both the solution (only the modes in the 
760 - 960 /iHz interval) for the whole dataset (continuous 
line) and the separate solution for the 4th subset (dotted 
line) are compared. The latter slightly increased the ampli- 
tude of the fit but in the maximum-amplitude cycle only 
about 50% of the amplitude is covered by the fit. Similar 
high- amplitude intervals are between BJD 3803.34 - 3803.44 
and 3850.46 - 3850.56. The inclusion of the rotationally split 
components helped to increase the amplitude but it was still 
not enough to match the observed amplitude. It is hard to 
imagine such a missing pulsation mode that could help to fit 
the high-amplitude cycles without destroying the acceptable 
fit of the low-amplitude cycles. It is more probable that some 
additional physical process is superimposed on the pulsation 
creating high- amplitude phases. Maybe we are faced with an 
extra effect of convection with the pulsation. 

A second characteristic feature is presented in panel b 
of Fig. [7] alternating high- and low-amplitude cycles fol- 
low each other. Neither the frequency solution of the whole 
dataset, nor the 4th subset can fit the light curve of these 
cycles. Similar features can be found between BJD 3847.51 
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iPfeiffer et all (|l996l ) determined Prot = 2.3 ± 0.3 d, 
also derived from an asymmetric triplet with {5f) — 
2.5 /xHz. This rotation p eriod agrees with our result within 
la. iHiirkal et alj (|2005l ) found ~3/iHz frequency splittings 
{{5f) — 3.27 ^Hz), corresponding to Prot ~ 1.8 d. Regarding 
that at each epoch different frequency content were deter- 
mined, the similar frequency spacings of the multiplets sug- 
gest they really correspond to stellar rotation at a constant 
rate. 

Beside the direct determinations, we searched for char- 
acteristic frequency spacing values applying a more sophis- 
ticated method. We selected the five main frequencies in 
the 65 - 85cd~^ range and - during a pre- whitening pro- 
cess - the highest amplitude ones in their vicinity. Then we 
performed the Fourier analysis of the frequencies obtained 
this way. That is we searched for regular spacing value(s) 
between the frequencies which may correspond to the rota- 
tional splitting phenomenon. To check our findings, we anal- 
ysed not only the whole frequency list, but some subsets of 
the frequencies, too. Our results suggest c haracteristic spac- 
ing values being around 2.6 and 3.7 /xHz (JBognar fc Papajol 
[2013) • These values are also close to the directly determined 
and the previously detected ones. 



Figure 7. Special features in the light curves of GD 154, a: in- 
tervals containing 4-6 cycles wliere the measurements represent 
much higher amplitudes for the cycles than the fit (maybe inter- 
action of the pulsation with convection), h: alternating high and 
low amplitude cycles follow each other (maybe chaotic behaviour 
of the pulsation). Fits for the solution for the whole dataset (con- 
tinuous line) and for the 4th subset (dotted line) are also given 



- 3847.56, 3851.48 - 3851.54 and 3935.34 - 3935.39. The 
presented case in panel b shows a phase shift between the 
solution for the whole dateset and the 4th subset. The al- 
ternating high and low amplitude cycles remind us of the 
chaotic behaviour of the stellar pulsation. 



3.3 Rotational splitting 

Thanks to regular observations over the whole season, we 
could find not only the normal modes but members of rota- 
tionally split frequencies. The direct determinations of the 
rotational splittings are 4.01, 2.99, 2.91 and 3.88 /iHz for 
the doublets presented in Table |3l respectively. Considering 
the three triplets listed in Table [5J we recognize that they 
are asymmetrically spaced and the tti = to -f-l splits are 
always smaller. This suggests that the « 920 /xHz doublet 
peaks, as also having small spacing, belong to m, = and 
m = -f 1 modes. 

We use the term 'triplets' for the multiplets we found, 
but without any mode identification, we cannot say which 
ones are real triplets (dipole modes) or three of five possible 
components of quadrupole modes. The asymmetric struc- 
ture of these closely spaced modes also makes difficult to 
distinguish between the I = 1 and / = 2 ones. Assuming 
that the triplet components at the dominant mode are high- 
overtone (A: ^ 1) Z = 1 ones, calculating from the 2.76 /xHz 
average split value, the rotation period of the star is 2.1 d, 
assuming solid body rotation. 



4 INVESTIGATION OF THE MAIN STELLAR 
PARAMETERS 

The efforts to determine the precise frequencies of the in- 
dependent modes were aimed at investigating the interior 
structure of the star. Asteroseismology gives us the oppor- 
tunity to provide constraints on the structure of the core, 
the hydrogen/helium layers, the mass of the star and to 
estimate the star's distance. For this purpose, we ran the 
White Dwarf Evolution Code (WDEC) originally written 
by Ma r tin Schwarzschild and mod i fied bvlKutter fc Savedof 




Lamb fc van HornI ||T97|V 'Winget' (119811). 'Kawalej 
WoodI (|1990D. iBradlevI |l9 93), Montgomery (1993^ 



and iBischoff-Kim. Montgomery fc Wingetl (|2Q0a ) . 

The WDEC evolves a hot (~100 000K) polytrope 
starter model down to the temperature we require, and 
gives an equilibrium, thermally relaxed solution to the stel- 
lar structure. For this model, the possible pulsation peri- 
ods of m = are determined by solv ing the non - radia l, 
adiabatic stella r pulsation equations IjUnno et al.l [l989|). 
iMetcalfd l|2001f) created an integrated form of the evolu- 
tion/pulsation codes, which allow us to obtain the period 
values with only one command. This way we can build model 
grids consisting of thousands of models in a very efficient 
way. 

We u sed the equation-of-state (EOS) tables 

of iLambl (Il974l') in the c ore, the EOS tables of 

ISaumon. Chabrier fc van HornI (|l995l ) in the envelope 



Iglesias fc Rogera 



Itoh et all (|l98l 



of the star, OPAL opacities updated by 
lll996l) and the conductive opacities by 
ll98J). The WDEC treats the c onvection by means of the 
mixing length theory (MLT) of iBohm fc Cassinellil (Il97ll ) 
using the a — 0.6 parametrization according to the model 
calculations of lBergeron et al.l (|l995l ) . The hydrogen/helium 
transition zone was treated by equilibrium diffusion cal- 
culations, while the helium/carbon transition layer was 
parametrized. 
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Table 4. Wc varied the stellar parameters and built model grids according to the minima and maxima values and step sizes given in the 
table. 



Grid 



Teff (K) 



M, 



{M&) 



-logAfn 



-logMn 



^O 



Xf, 



1 


10 600-11800 


0.600-0.800 


2 


4-11 


2 


10 600-11800 


0.600-0.800 


2-3.5 


4-11 


Step sizes: 


200 


0.005 


0.5 


0.2 



0.5-0.9 0.1-0.5 

core profiles by ^alaris_e^^L 

0.1 0.1 



We fitted the observed pulsation periods with the cal- 
culated dipole and quadrupole ones and searched for the 
best-fitting models. The goodness of the fit was described 
by the r.m.s. value calculated by the following way: 



EL(^^ 



iV 



(1) 



where N is the number of observed periods. Th e ar.m.s. va l- 
ues were calculated using the fitper program of |Kim| (|2007J) . 

We built two model grids using different core composi- 
tion profiles. In one case, we varied five input parameters of 
the WDEC: Te«, M,, Mh, Xq (central oxygen abundance) 
and Xfm (the fractional mass point where the oxygen abun- 
dance starts dropping) and fixed the mass of the helium 
layer at the 10~^ M* 'canonical' value. In the second scan, 
we varied the mass of t he he lium layer, but used the core 
profiles of ISalaris et al.l l| 19971 ) based on evolutionary calcu- 
lations. For this second grid, only four stellar parameters 
were scanned. Table |4] shows the parameter space covered 
by our grids and the step sizes applied. 

The Tcff and logg values of GD 154 determined by 
high signal-to-noise optical spectrophotometry a re 11 180 K 
and 8.15 dex, respectively (jBergeron et al.lll995^ . Consider- 
ing that the external uncertainties are estimated to be hun- 
dreds of Kelvins and could achieve ±0.1 dex, we decided to 
cover a relatively large range in effective temperature and 
surface gravity. A DA white dwarf with logg = 8.15 has a 
mass of ~ 0.7 Mq. In our grids, we searched for the best- 
fitting mod els between log.*? ^ 8.0-8.3 (0.6-0 .8 M(7^, see 
the tables of iBradlevI 19961) . iBergeron et al.l (|2004l ) estimated 
~ 200 K and ~ 0.05 dex for the external errors of Tch and 
log <;, respectively. According to these values, our grid covers 
±3(T range in Toff and M«. 

Pfciffer et al. ( 1993) suggested a very thin hydrogen 
layer for GD 154 with the mass of 2(±1) x 10~"' M,. Tak- 
ing their result into account, we allowed for hydrogen layer 
masses between 10"* and 10~^^ M, in the grids. 

We also used the database of ZZ Ceti periods (dipole 
and quadrupole) calc ulated from fully evolutionary models 
(jRomero et al.ll2012al ). The period values were derived by 
models with consistent chemical profiles from the core to the 
surface. The authors allowed the evolution of the stars from 
the zero-age main sequence in the calculations of these pro- 
files. More details on the code, the input physics and a num- 
ber of exa mples of its asterosei smol ogical applicat i ons ca n 
be found in lAlthaus et al.l (|201G| ) and lRomero et al.l(|2012bl ). 



4.1 Parameters of the best- fitting models 

We used two slightly different sets of periods for the as- 
teroseismological investigations of the star. Our frequency 
analyses and tests show that the periods given in Table [3] 



describe the light variations of GD 154 in the 2006 observa- 
tional season well. Accordingly, the values used in our first 
run are: 402.6, 1088.5, 1131.8, 1160.7, 1191.7 and 1238.2s, 
respectively. The second run differs from this in one period 
only: we replaced the 1191.7 s mode with 1186 s, the one 
with the second largest amplitude. During the 1991 WET 
observations, the period of the dominant mode was 1186.5 s 
and was found to be the m = component of an (asym- 
metric) triplet (jPfeiffer et al.l ll996l). As Tables and E] show 
in our dataset the 1191.7 s period dominates. This one is 
still part of a triplet structure and could be regarded as 
the m — —1 peak of the 1186 s mode. We used two differ- 
ent sets of periods because of this ambiguity. This way, we 
could examine the effect of slightly different periods on the 
best-fitting models' parameters. 

Even though we doubled the number of known period 
values applied in the star's seismic investigations compared 
to previous research, GD 154 is still not a pulsator rich 
in known modes. This means that we find several models 
with low ar.m.s. values as the result of the fitting procedure. 
Therefore, we applied further constraints during the model 
selection: assuming better visibility of Z = 1 modes, we se- 
lected the models with at least three I — 1 solutions to the 
observed periods. As an additional criterion, we considered 
the ones which give I — 1 value for the dominant (1186 or 
1191.7 s) mode. However, as the average period spacing of 
the long-period modes is low (below 40 s), we assumed that 
the observed modes are not solely I — 1 ones, because this 
would require a much higher stellar mass for GD 154 than 
the spectroscopic value. 



4-. 1.1 Stellar parameters 

The left panel of Fig. [8] shows the 6-1-6 best-fitting mod- 
els in the Tcs-M, plane for both period lists. The num- 
bers indicate the models' hydrogen layer ma sses and we also 
denot ed the spectroscopic solution given bv I Bergeron et al.l 
([1993) with it's uncertainties. As it can be seen, we find the 
best-matching (three) models for the spectroscopic mass us- 
ing the second period list. Table [5] (rows 1-3) summarize 
the parameters of these models. In two cases (M* — 0.71 
and 0.73 Mq) the mass of the hydrogen layer is around 
lO'^'M*, while the model star with 0.705 Mq has Mh = 
1.6 * 10~*M, , which means a considerably thinner layer. 
The best-matching model to the spectroscopic parameters 
is the one with M, = 0.73 Mq and T^h = 11 200 K. In this 
case the mass of the hydrogen layer is Mh ~ 6.3 * 10~^ Af, . 
When we investigate a pulsator showing amplitude vari- 
ations, it is worth checking if we can add further modes to 
our period lists observed i n a different s eason . We found 
one mode, also presented bv lHiirkal et al.l (|2005l ): the one at 
1271.5 s. Then we have four period lists: two with 6-6 and 
two with 7-7 periods. In the latter case we selected the mod- 
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Figure 8. Plots of the models that best fit the data for six (left panel) or seven (right panel) modes in the T^ff — M* plane. Circles and 
squares denote the solutions obtained by including the 1191.7 or 1186 s mode into the period list. Their ar.m.s. values are colour-coded. 
The models' hydrogen layer masses (-logAfjj), the spectroscopic T^g values, and their uncertainties are also indicated. The background 
grid corresponds to our model grid's step sizes. 

Table 5. Parameters of the selected models using different core profiles and slightly different period values. The observed T^g and J\/» 
values IIBereeron et al.lll995r) are are given in the last row. 



No. Tefl 



M* 



-logA/Ho -logMn Xo Xf, 





(K) 


(Mo) 




1 


11600 


0.710 


2.0 


2 


10 800 


0.705 


2.0 


3 


11200 


0.730 


2.0 


4 


11000 


0.720 


2.0 


5 


11600 


0.710 


2.0 


6 


11200 


0.675 


2.0 


7 


11400 


0.675 


2.5 


8 


11000 


0.710 


2.0 


9 


11200 


0.700 


2.0 


10 


11000 


0.710 


2.0 


11 


11241 


0.705 




12 


11639 


0.705 




Observations 








11180 


0.70 





Period values in 

{l,k) 



seconds 



(s) 



6.0 



6.2 



5.8 



6.0 



5.2 



5.2 



5.6 



4.2 



5.6 



4.445 



9.339 



0.8 



0.6 



0.6 



0.8 



0.2 



0.4 



0.4 



0.8 0.5 



0.2 



0.73 0.54 



0.73 0.54 



0.73 0.54 



0.73 0.54 



0.73 0.54 



402.0 


1088.9 


1132.6 


1159.2 


1184.9 


1238.4 




(1,7) 


(1,23) 


(2,42) 


(2,43) 


(1,25) 


(2,46) 




402.9 


1087.6 


1131.5 


1162.0 


1188.1 


1238.3 




(2,11) 


(1,19) 


(2,35) 


(2,36) 


(1,21) 


(1,22) 




402.0 


1090.8 


1131.2 


1161.3 


1185.8 


1239.2 




(1,6) 


(1,22) 


(2,40) 


(2,41) 


(1,24) 


(2,44) 




400.4 


1088.5 


1132.2 


1159.5 


1194.2 


1238.2 


1273.4 


(1,6) 


(1,22) 


(2,40) 


(2,41) 


(1,24) 


(1,25) 


(2,45) 


402.0 


1088.9 


1132.6 


1159.2 


1184.9 


1238.4 


1268.5 


(1,7) 


(1,23) 


(2,42) 


(2,43) 


(1,25) 


(2,46) 


(1,27) 


402.4 


1089.2 


1133.6 


1158.1 


1191.9 


1240.1 




(2,12) 


(1,22) 


(2,41) 


(2,42) 


(1,24) 


(1,25) 




402.2 


1086.5 


1132.7 


1160.5 


1194.5 


1239.4 




(1,6) 


(1,22) 


(2,41) 


(2,42) 


(1,24) 


(2,45) 




398.9 


1085.6 


1131.4 


1160.7 


1188.0 


1236.8 




(1,6) 


(1,22) 


(2,41) 


(2,42) 


(1,24) 


(1,25) 




398.9 


1085.9 


1132.8 


1159.9 


1192.7 


1236.0 


1269.9 


(1,7) 


(2,43) 


(2,45) 


(2,46) 


(1,27) 


(1,28) 


(1,29) 


398.9 


1085.6 


1131.4 


1160.7 


1188.0 


1236.8 


1267.7 


(1,6) 


(1,22) 


(2,41) 


(2,42) 


(1,24) 


(1,25) 


(2,46) 


399.3 


1082.6 


1133.9 


1160.0 


1192.3 


1239.8 




(1,7) 


(2,43) 


(2,45) 


(2,46) 


(1,27) 


(1,28) 




404.3 


1085.3 


1130.7 


1165.7 


1182.6 


1237.9 




(1,5) 


(1,20) 


(1,21) 


(2,38) 


(1,22) 


(1,23) 




402.6 


1088.5 


1131.8 


1160.7 


1191.7 
1186.0 


1238.2 


1271.5 



0.86 



1.09 



1.13 



1.51 



1.38 



1.52 



1.55 



2.16 



2.07 



2.45 



2.97 



2.92 



els with at least four I — 1 modes instead of three. As the 
right panel of Fig. [8] shows, the parameter space occupied 
by the best fitting seven-period models differ only slightly 
from the six-period ones. Using seven periods most of the 
models have M* > 0.7 Mq. Considering the hydrogen layer 



masses, the average values are ~ 10"'' and 4 * 10"'' M* for 
the six- and seven-period solutions, respectively. We can rec- 
ognize in Fig.[8]that there are common models, the ones with 
Toff = 10 600 and 11 600 K and Mh = 10"^" and 10"*' M*. 
Our best-matching models using seven periods are the 0.720 
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Figure 9. Models in the Mt-Mjj plane with (Jr.m.s. < 2.5s 
and Salaris et al.'s core profiles. The six-period light curve solu- 
tion with the 1186 s mode was used. The models' rrr.m.s. values 
are colour-coded. The background grid corresponds to our model 
grid's step sizes. The figure shows the trend that more massive 
stellar models have thinner hydrogen layer. 



and 0.710 Mq ones with Tcft = 11 000 K and 11 600 K. They 
have Mh = 1.6 * lO"'' and lO"'' M, (Table [3 rows 4-5). As 
can be seen, models No. 1 and 5 have the same parameters, 
this is one of the common points. 

Considering the solutions with Salaris et al.'s core pro- 
files and within the la limit in M*, the hydrogen layer 
masses are between 10~* and 10~® A/* in most cases and 
the Mhc = 10"^ M, value is preferred. Rows 6-10 of Table[S] 
show some of their parameters; they are selected on the basis 
of having T^s close to the spectroscopic value. These models 
have Mh = 6.3 * 10~^ 2.5 * 10"^ or 6.3 * 10"^ M*. 

Investigating the parameters of the best-fitting models, 
we found a trend that more massive stellar models have a 
thinner hydrogen layer. An example for the phenomenon can 
be seen in Fig. [51 where we plotted the best-fitting models 
with Ur.m.a. < 2.5 s iu the M, - Mh plane. The explanation 
of this phenomenon is that generally for thinner hydrogen 
layers, the mode trapping cycle is longer and the average pe- 
riod spacings between the consecutive overtone modes are 
also larger (i.e. there are larger differences between period 
spacing minima and maxima). This can be partially offset 
by the shorter average period spacing of a higher mass star 
l|Bradlevl[l996t ). The panels of Fig. [TD] demonstrate the in- 
fiuence of the varying stellar parameters on the period spec- 
trum of a model star. The stellar mass and the mass of the 
hydrogen layer has a great effect on the observed periods. 
We can find a similar period structure for the lower stellar 
mass - higher hydrogen layer mass and higher stellar mass 
- thinner hydrogen layer models, and this results in the ob- 
served trend seen in Fig. [9] 

Using the coarse grid of iRomero et al.l (|2012al ) we ob- 
tained the best-matching model to spectroscopy applying 
the first period set. It has Toff = 11 241 K, M, = 0.705 Mq 
and Mh = 3.6 * 10"^ M, (Table [3 row 11). The second 
best-matching model belongs to the second period list and 



has higher effective temperature (11 639 K) but the same 
mass (Table [5] row 12). The largest difference is in the 
mass of the hydrogen layer, which is significantly lower, 
only 4.6 * 10"^'' A/*. Considering all the solutions within the 
lyr.m.s. < 3 limit, the mass of the hydrogen layer is between 
3.5 and 7.6 * 10~^ M, , from which the latter one is the only 
exception. 

Knowing the luminosities of the selected models, we es- 
timated the star's asterosei smological dis tance by calculat- 
ing the distance modulus (|Bradlevll200ll ). The log(I//I/Q) 
values of the models in Table [5] are between —2.7 and —2.8. 
The seismological distances calculated from the ten differ- 
ent models of Table [S] are between 41.5 and 46.4 pc. The 
average distance and parallax value is 44.2 pc and 22.7 mas, 
respectively. This result is close to the 45.4 pc value derive d 
by spectroscopic observations of iLaioie fc BergeronI (|2007|), 
which shows that the pulsation analysis supports the spec- 
troscopic result on the distance parameter. 

For the sake of completeness, we summarized the results 
obtained by previous asteroseismological investigations of 
GD 154 in Tabled 



4-1.2 Mode identification 

As it was already mentioned, we placed a constraint only 
on the I value of the dominant mode: we assumed that it is 
I — 1. All of the other modes were allowed to be ? = 1 or 
2. Considering the I and k values of the selected models in 
TableO we obtained the same I only in the case of the 1160 s 
mode. This result suggests that it may he an I = 2. None 
of the selected models gives I = 1 solutions for all of the 
periods. We cannot uniquely assign an I value for the other 
modes, however, the 402, 1088 and 1131s modes' I values 
are 1, 1 and 2, respectively, in the vast majority of cases. 
Assuming that the 1131 and 1160 s modes are I = 2, they 

represent consecutive overtones. 

The models selected by ICast anhcira fc Keplen (|2009l ) 
and lRomero et al.l (|2012bl ) (see Table [6]) give I = 1 solutions 
for the 402, 1088 and 1186 s modes as weU. 



4.1.3 Mode trapping and the mass of the hydrogen layer 

IPfeiffer et al.l ([1996) assumed that mode trapping - as an ef- 
ficient mode selection mechanism - could be responsible for 
the small number of observed modes in GD 154 and deter- 
mined a very low mass for the hydrogen layer (~ 10~^° M,). 
As Table [S] shows, the results of the previous asteroseismo- 
logical investigations also supported the presence of a very 
thin hydrogen layer. However, considering our selected mod- 
els in Table \S\ we found that the masses of hydrogen lay- 
ers were higher than expected. Our model solutions have 
Mh values between 6.3 * 10~^ and 6.3 * 10~^A/* (with an 
average value of 3.9 * 10~®) except for two cases, when 
Mh = 1.6 * 10"* and 4.6 * 10"^° M*. However, these lat- 
ter models are below or above the la limit in effective 
temperature. Our results show that although GD 154 may 
have thinner hydrogen layer than the maximum allowed 
(A/h ~ 10~*), no extremely thin layer is necessary to ex- 
plain the sequence of the observed periods. 

The question is raised: why the other authors obtained 
lower mass values for the hydrogen layer. The answer may 
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Table 6. Physical parameters and mode-identification results pr esented by diff erent authors. They selected the model solutions using 
the period values determeined by the 1991 WET campaign (,Pfciffer et al.lll99(i ). The atmospheric parameters and periods obtained by 
observations are denoted by asterisks. 



Teff (K) 


M, (Mq) 


-log Mho 


-logAfn 


Period values in seconds 


Rcf. 


11 180* 


0.7* 




9.7 


402.6*{1) 


1088.6* (2) 


1186.5* (1) 


Pfeiffer et al. (1996) 


11200 


0.68 


2.0 


7.5 
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1088.5(1,19) 


1186.9(1,21) 


Castanheira & Kepler ('2009) 


10 800 


0.73 
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Figure 10. Period spacing diagrams of model stars with T^g = 11 200 K, Afjjc = 10~^ A/* and 50/50 C/O fraction core. Filled and open 
circles denote the I = 1 and 2 modes, respectively. The panels show the influence of the stellar and hydrogen layer mass variations on 
the period spectrum. 



be the number of m o des ii sed in the period fi t s. Both 
ICastanheira fc Kepler! l|2009l 'l and iRomero et al] ll2012bl) 
worke d with the three periods determined bv lPfeiffer et al.l 
l|l996t ). This fact obviou sly strongly in f luenc ed the outcome 
of the fit. The results of iPfeifFer et al.l l|l99a ) have not been 
obtained by performing period fits with a model grid. Be- 
cause of the small number of known modes they needed to 
impose further constraints. As mode trapping means an effi- 
cient mode selection mechanism, they assumed that the ob- 
served modes are trapped in the outer hydrogen layer and 
derived the Mh ~ 10"^'' M, value. In our case, the increased 
number of observed modes is sufficient for the asteroseis- 
mic investigations using model grids, including searching for 
trapped modes among the observed ones. 



We constructed the (forward) period spacing diagrams 
for our selected models (Fig. Ill[) . Minima in the diagrams 
denote departures from th e uniform period sp acings caused 
by mode trapping (see e.g. lBradlev fc Wingetlll99lh . Fig.fTTl 
shows that the long-period modes do not show tendency 
to occur around minima, so our model selection does not 
support the idea of mode trapping at long periods as a mode 
selection mechanism. The best candidate for a trapped mode 
is the 1238 s one according to Fig. 1111 However, the situation 
is different in the case of the only short-period mode. In most 
cases we can find it near or at a period-spacing minimum. 
This suggests that the 402 s mode could be a trapped mode. 



5 SUMMARY AND CONCLUSIONS 

Our comparative period search confirmed the frequencies 
found previously, except the one at 787 //Hz given by 
iHiirkal et al.l (|2005l ). Additionally, we localized modes at 
807.62 and 861.56 /xHz values that have never been reported 
before. We confirmed by test investigations, that six modes 
can be considered as independent normal modes of pulsa- 
tion. 

Four doublets around the largest amplitude modes were 
directly found. The dominant mode of our whole data set 
(839.14 /iHz) and the dominant mode observed in two previ- 
ous pulsation stages (843.15 /iHz) are members of the same 
rotational triplet. This latest member of the triplet became 
dominant at the end of our observing run, when GD 154 
also presented a monoperiodic pulsational stage. 

We localized the second harmonic of the mother fre- 
quency at 2517.48 /tHz, near the high frequency normal 
mode at 2484.14 /iHz, and some linear combinations too. 
However, no subharmonics reported by previous observa- 
tions were found over our whole observational season. Char- 
acteristic features are localized in the light curves partly 
suggesting an effect of convection to the pulsation and re- 
minding us of the chaotic behaviour of stellar pulsation. 

Comparative analyses of subsets revealed and test inves- 
tigations confirmed a remarkable intrinsic amplitude change 
of frequencies at 839.14 and 861.56 /iHz, although part of it 
can be caused by unresolved rotational triplets. 
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Figure 11. Period spacing diagrams of the ten different mod- 
els of Table [5] Filled and open circles denote the Z = 1 and 2 
modes, respectively. Vertical dashed lines mark the observed pe- 
riod values. We denoted the periods analogous to the observed 
ones in the given model with open squares. The models' numbers 
corresponding to Table [5] are also indicated in the panels. 



With our new frequencies we have doubled the number 
of period values applied in the star's previous seismic inves- 
tigations, which allows a more detailed study of GD 154. We 
found models with effective temperatures and masses within 
the 1(7 limit of the spectroscopic values (« 11 000 — 11 400 K 
and 0.68 — 0.73 M*) that fit the observed periods well and 
also give / = 1 solutions for at least half of the modes. 
The best-fitting models have hydrogen layer masses between 
6.3 * 10~^ and 6.3 * 10~^ M*, which suggests orders of mag- 
nitudes thicker layer than previously published. The expla- 
nation of this difference may be the number of modes used 
for the seismic studies and that we did not assume that our 
observed modes were trapped ones. Considering the mass 
of the helium layer, our results also supports the 10~^ A/hc 
'canonical' value. 

The known luminosities of our selected models allowed 
us to determine the seismic distance of GD 154. In agree- 
ment with other authors' finding, the average value calcu- 
lated by our selected models is 44 pc. We also investigated 
the possibility of mode trapping, constructing period spac- 
ing diagrams. Our results do not support the idea of mode 



trapping at long periods as a mode selection mechanism, and 
suggest that the short-period mode may be a real trapped 
mode. This result shows that we do not have to presume 
that all the observed modes are trapped ones. 

Both the new frequency content, the altering pulsation 
between the multi- and monoperiodic stage and also the 
need for more constraints for modeling express a require- 
ment for further investigation of the complex behaviour of 
GD 154. Regular monitoring of the star would be necessary 
not only to detect new pulsation modes for asteroseismic 
investigations but to follow-up more mono- and multiperi- 
odic stages. This might allow us to determine whether this 
altering between the two pulsation stages has any regular 
behaviour. 
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